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ABSTRACT 

A mathematical model is developed for determining an optimum 
submarine tender load list for the POLARIS (FBM) weapon system. 
Military essentiality of spare parts and the stowage capacity of the 
tender play prominent roles in the development of the model. The 
measure of effectiveness which is maximized is the average number 
of demands filled over a given period of time. This is equivalent to 
minimizing the number of stockouts of an item. An example of the 
use of the model is included for the case where the demand for 
items is characterized by the normal distribution. The paper in- 
cludes a computer flow chart and program to aid in implementing 
the model. 

The authors wish to express their appreciation for the 
assistance and encouragement given them by Associate Professor 
Julius H. Gandelman of the U. S. Naval Postgraduate School. 

The authors are also indebted to Miss Patricia Hoang for her assis- 


tance in programming the model on the CDC 1604 computer. 
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GLOSSARY OF TERMS 


a’ given item where i=1,2,3...... .N (the entire 
list of spare parts is drawn up and then a number 
is assigned to each spare part to identify it). 


amount of item ito be stocked. 


-th 


demand for the i item. 


probability that an amount x; will be demanded 
during the period under consideration. 


expected demand for item i ; = mj. 


expected shortage of item i. 


sth 


average number of i” item supplied to meet demand. 


Lagragian multiplier. 


th item. 


cube (volume) of the i 
total cube available for spare parts stowage on tender. 


total cube associated with any X for which = CLZ.FCL. 


essentiality (military worth) number assigned to a 
item. 


ee aL : 
deviation of i - item from the mean, 


mean of the qe item. 


characteristic number associated with the jth item. 








1. Introduction. 

The trend toward the development and utilization of highly com- 
plex, sophisticated weapon systems within the Navy has generated 
an ever-increasing requirement for more efficient logistical support. 
This paper is an attempt to develop a mathematical model that can 
be used as a decision aid for one aspect of the logistics problem 
associated with the POLARIS system. The model provides a method 
for determining the items and associated quantities that should be 
stocked aboard the supporting tender. 

The requirement for more effective tender load lists is even 
more evident when viewed in conjunction with the POLARIS mobile 
base concept. 

In the event of enemy attack, the Navy's immediate POLARIS 
retaliatory capability will have been dependent on the ability of the 
supporting tender to provide the necessary SSB(N) replenishment 
and repair requirements. This aspect of the tender function is fur- 
ther complicated by the limitations of stowage space available for 
stocking the required spare parts. Asa result of this space 
limitation, tender load lists must wetlect SSB(N) requirements, and 
tender re-supply procedures must be closely coordinated and 
evaluated to insure that SSB(N) requirements are being met in an 


effective and efficient manner. 
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The concept of military worth is incorporated in the model to 
reflect the assignment of priorities between the various items. 
These essentiality numbers are an attempt to define the military 
worth of a particular item with reference to the performance of the 
integrated weapon system. Essentiality factors will be assigned 
arbitrarily; however, it should be apparent that actual MEC numbers 
could have been utilized for illustrative purposes. 

The measure of effectiveness to be maximized will be the average 
number of demands for spare parts that are filled by the tender over 
some period, as a function of the military essentiality of the various 
items. This is equivalent to minimizing the number of stockouts 
Of an item. 

It is obvious that if the demand for an item was known exactly, 
the problem would be simplified considerably subject to the space 
limitation of the tender. Although demand cannot be predicted with 
certainty it may be possible to estimate the probability distribution 
of demand for each item on the basis of past usage data. In the event 
that usage data are not available, the model can be utilized by as- 
suming a probability distribution of demand for the item(s), based 
upon experience and judgment. This would be the situation encoun- 
tered in developing an initial tender load list for a new class of 


submarines. As usage data become available the tender load list 
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can be continually updated to reflect actual demand, utilizing 
automatic data processing equipment. Appendix III contains a com- 
puter flow diagram and program that will accomplish this objective. 
Although this paper is primarily concerned with presenting a 
method for determining an optimal tender load list, certain aspects 


of the tender re-supply problem are included for consideration. 





ap Assumptions 

The following simplifying assumptions are introduced in order 
to illustrate the proposed method. 

1) The model is concerned with spare parts and does not 
consider dry stores and other consumable items. Asa result, C 
refers to the cube available for spare parts stowage. 

2) x; 18 a continuous variable, and its associated probability 
density function f;(x;) possesses the properties required for the 
existence of derivatives and integrals. 

3) The demand for each item is independent of the demand 
for every other item. 

4) The probability distribution of demand is known for each 
item over the time period considered. 

5) The list of items considered reflects only those for which 
demand has been generated. If demand for an item stops, then the 
item is removed from stock and a new load list is generated. Con- 
versely, if demand for a new item is generated, it is added to the 
load list and a revised Z=|z, Zoy+-ZN is determined by means 
of the model. 

6) The development of the model is independent of the re- 
supply phase of the tender operation. The load list determined is a 
function of demand over some period during which the tender must 


provide support without being re-supplied. Theoretically the tender 





is completely depleted of spare parts at the end of the period con- 
sidered. The effect of re-supply as a function of the desired supply 
objective will be considered in a later section. 


7) The utility values defined in Appendix I are additive. 








oO. Presentation of the Model. 

This section presents the essential concepts of the proposed 
model. Appendix I contains a detailed mathematical formulation of 
the results of this section. 

The expectation, EU J], of x. with respect to the 
probability function P+] gives the following expression for the 
average amount of item i demanded: 


oo 
(1) Bex = aa x; £; (x,) dx;, 


h item. 


where x;= demand for the it 
(The limits of integration range from Q to©O since only items 
for which a demand is generated are considered) 

Consider the following situations: 

(A) x <3 s 
(b) Xi > Zj: 
where Z; = amount of item i to be stocked. 

If demand for an item is less than the amount of the item that 
isstocked then the tender can meet demand. However, if demand 
exceeds the stock level then only an amount zj can be supplied. If 
this condition occurs (x; - 2;) represents the shortage (backorder) 


associated with that item. The expected shortage of an item can be 


expressed as a function of the amount of the item stocked. 


C8 
(2) bi Cz = if (xj— 23) £,(xj) dx; For x;> 2Zj, 
£5) 
= 0 otherwise . 





The expected value of the amount actually supplied to meet 


demand is given as a function of the amount stocked on the tender. 
This results in the following expression: 
ai co 
(3) Sylec ee { , x; f;(x,) dx; + z, he. £; (x; )dx; 
Manipulation of equations (1-3) confirms the intuitive concept 


that the expected demand must equal the expected demands filled 


plus the expected shortages. i 


(4) Ex = siedt+bi (29 


As mentioned previously, the measure of effectiveness is a 
function of the military essentiality of the item as well as the de- 
mands filled by the tender. This requirement is met by defining 
M(z):Z5; 2+ Zy)= 2M, S;(z;) in accordance with assumption (7). 
The problem is now one of determining the set Z =|] Z,,Z2,+--Zeal 
such that M UZ ee ee) = = mM Seay } 


is maximized subject to the constraint = c: zi =C, ZnO 
GA 


Equation (5)(below) is derived in Appendix I and is a general 
result of the considerations mentioned in this section. This equation 


provides the means for determining an optimum tender load list. 


z: 
(5) J fil ag = 1- ACER), ie 2. 


L 
where A is a lagrangian multiplier with 


values inthe range min ( ws = N2O 





CL ~— 
subject to: NI 
2 auz=C = Liane 


Appendix I 





e «(eee «<4 





> eet" 
- —— te 
. 4 ee ee 
9 ae 


_—— 


4. General procedure for utilizing the model. 

The general procedure for solving equation (5) is to assume a 

value for : in ( Ss 

ue for corresponding to min ( <=), and compute the set Z 
by using the appropriate cumulative distribution function table. The 
C' obtained (where C' results from using a A other than that for 

WN 

which 2.¢i2.=C) when the z; are substituted into the constraint 
NI ' 

equation ( a C.Z.= C ) determines the minimum required stowage 
Space aboard the tender as a function of demand over the demand period 
associated with f;(x;). (As X increases, C' decreases. Since 

Mia... sie : 
min ( q+) is the largest permissible value of A , the corresponding 

Nee 

C/is a minimum.) It should be noted that for any A in the range 
min Se AZO , every 2; 27 0. This automatically sets a mini- 
mum stowage requirement on the tender of at least the sum of the 
cube of each item. 

Under the assumption that C is known for the tender, a com- 
parison is made between this value and C'. If C > C' an iterative 
solution is required to determine the X which results in C'=C. 

The corresponding z;'s obtained are the optimum load list. 

If, on the other hand, C < C' the model indicates that the ten- 
der cannot effectively meet demands without re-supply during the 
period under consideration. In the event this occurs several decision 
criteria can be adopted. If demand for each item is relatively uni- 


form, one solution would consist of decreasing the tender load list 
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time requirement. > For example: Let C' = 15,000 cu.ft.represent 

a ninety-day load list minimum requirement. If C = 5,000 cu. ft., 
then the tender can in effect maintain a 30 day load list by decreasing 
each 2; proportionally. | 

This criterion has the advantage of maintaining a load list which 
consists of all items that are demanded. The most obvious disadvan- 
tage is the increased re-supply requirement. 

Another possibility is the elimination of items from the load 
list. This requires a decision criterion to select the item(s) which 
should be deleted. Since C' decreasesas XA increases.the model 
indicates that the item corresponding to min( Bs ) should be the first 
item Resear from stock. If (C'—c, a, =C; a. < Z.) then an 
amount (z; -a;) can be maintained for this item. If(C'—cj2;>C) 

Mie 


then the item corresponding to min (-<* ) is eliminated completely. 


A new min( a ) is selected from the remaining list of N-litems 





and the process is repeated. Whena A can be found such that C=C! 
the z's have been determined for this decision criterion. This 


particular policy insures that the items with the highest essentiality 


Mc 


' Z 
to cube ratio ( = ) are retained on board. 





That period during which the tender must be able to operate 
without benefit of re-supply and still maintain effective support 
capability. 


¢ Note that other policies could have been chosen. For example, 
one might choose a policy that insures that items with high values of 
P are retained on board, while items with low values of P are eliminated. 





It permits tender support to be maintained at some reduced effective- 
ness level, (This reduced level can be a significant one), In effect, 
this is analogous to an unexpected demand situation and will be con- 
sidered further in connection with tender re-supply requirements. 
The iterative nature of the solution and the number of items 
comprising load lists indicates the necessity for utilizing electronic 
computers if the system is to be incorporated effectively. Appendix 


III contains a computer description and flow chart for the model. 
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ae Example and General Results: 

A numerical example will serve to illustrate the procedure 
for determining an optimal tender load list utilizing the proposed 
model. Consider the case where the demand for each item is 


characterized by the normal distribution. From equation (5) 


Zo Mc 


= ( en “Sd = BCBS) LAG), 
© 





\ar 


Subject to the available stowage constraint. It should be noted that 
each item could have different demand distributions and the model 
would still be valid. The only requirement for utilizing the model 
is that the density functions be continuous. 

Table (1) contains the parameter values pertinent to the 

example. In accordance with the procedures of the previous section, 
an available tender capacity of 10,152 is assumed. It should be 
apparent that due to the iterative nature of the solution, ay tte: 
will be approximately equal to 10,152 for an optimum set Z. 
(The most obvious reason is the need for rounding off the values 
obtained to intergers. In order to insure that the space constraint 
is met, the z;'s obtained should be rounded off to the lower interger 
with the result that - ez © y 

The values assigned to the items listed in Table (1) reflect 


varying combinations of essentiality, cube, expected demand and 


laAnother model could be developed which would utilize discrete 
functions. 


1] 
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standard deviation. 
The essentiality valves in Table 1 were arbitrarily chosen to 
show that items 1-8 are approximately four times as essential as 
items 9-16. Ina real life problem involving 5 large number of 
items, a scheme using the military essentiality codes could be adapted 


for use in this model. 
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In order to determine the value of M once an optimum set Z 
is obtained, equation (3) must be utilized. For the case when 


demand is characterized by the normal distribution, 


Zi oo 
(3) 3;C2.] = j xf. (xX: Jax. a Pils $x) Ax and 
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Appendix II contains a detailed derivation of equation (3'.). 
Since equation (5) requires z, > 0 this implies that for the 
case of the normal distribution equation (5) and (3') are only valid 
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> 3.87. (Fig. 1) 
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Table 1 indicates that this condition has been satisfied. 


Due to the requirement above, the third term of equation (3') 
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- The value fora is dependent on the accuracy of the tables 
used. The value of 3.87 was obtained from the Chemical Rubber 
Company Standard Mathematical Tables. 
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can be ignored in this example since its contribution is insignificant 
for the values assumed. 

Table 2is a presentation of the results obtained when the 
parameters shown in Table 1 are used inthe model. The optimal 
set Zo =|| ny tes. + + Ane || which most nearly satisfies the arbitrary 
constraint of 10,152 cu. ft. tender capacity was found by using A=,04\6. 

Note that items 9-12 generally lag in percent of expected average 
demand filled S:2Vekd, for each value of A. This is intuitively 
appealing because these are items with low military essentiality and 
high cube. The items with fen essentiality and low cube (5-8) are 
stocked at consistently high levels. 

Figures 2a thru 2d are plots of percent of expected average 


number of demands filled, versus a characteristic number P, defined 


as P, = Bee 
. oC 


Since each item is characterized by the parameters m,W, Mandc, 
each item is also characterized by P. Note that for different values 
of C' the curve shifts somewhat but retains the same general shape. 
In each case as the characteristic number increases, the percent of 
expected average number of demands filled increases. 

The effect of variability of demand for this example is investi- 
gated in Figure 3 by observing the behavior of items 9-12. Since 
items 9 and 10 (also 11 and 12) differ only in their values of 


we would expect any differences in their optimal stock levels (2Z() 
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for various total cube requirements to be due to variability of demand 
Ge) . 

Figure 3 supports the argument that when a large total cube 
is available, items with high values of SW will be stocked in larger 
amounts than items with small GS . However, when the total cube 
available becomes more restricted, the stock of items with high 
variability of demand is reduced rapidly. Plots for other items which 
differ only in their value of S show a similiar trend. 

Reference to Figure 4 indicates that as total available cube in- 
creases, the total essentiality of the load list increases also, but at 
a decreasing rate. Planners who wish to increase the spare parts 
capacity of a tender would have to weigh the cost of increasing 
capacity against the value of the added cargo. Use of a computerized 
version of the model with a large number of items would provide a 
new optimal load list for the proposed increased capacity. Reference 
to a curve similiar to that shown in Figure 4 would be an aid to the 
planner who must determine whether the essentiality of the added 


inventory is worth the cost of increasing tender spare parts capacity. 
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Fig 3: Effect of Variability of Demand 
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6. POLARIS Tender Re-supply (General) 

The maximum deterrent value of the FBM weapon system can 
be realized only if the maximum number of SSB(N)'s are maintained 
in the designated areas. The logistic support concept is based on 
satisfying this requirement. Lack of essential spare parts can be 
interpreted as a deficiency in the number of ready missiles scheduled 
to be on station. 

A load list is developed to allow the tender to provide effective 
support for a given period of time in the event re-supply is cut off 
for any reason. Asa result, replenishment of the tender is in- 
fluenced by the following factors: 

(1) maintain stock levels at or near the optimum load list 

point, and 

(2) provide delivery of nonstandard demands (items not 

included on the load list for one reason or another) 
from the stock point to the tender. 

The factors considered in selecting the mode of transportation 
(air or surface) for an item include the transit time available, 
priority, physical characteristics of the material, capability and 
availability of the various carriers, security and cost. Of these 
factors those necessary to effect delivery to meet program require- 
ments should govern and considerations of cost and equitable traffic 


distribution should be regarded as secondary. 
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In order to illustrate the considerations mentioned above con- 
sider the case where a tender is assigned to support X SSB(N)'s. 
The load list for the tender is designed to provide Y days support 
during which time each SSB(N) is alongside for 42 days. (Fig 5). 

Y can be interpreted as one cycle during which an SSB(N) has an 
upkeep period of % days and a patrol of W days (Y = #+W). If the 


assumption is made that the alongside availability period for each 





SSB(N) is the same, then A= =K is the number of SSB(N)'s at 
the tender at any time. Under this as sumption an SSB(N) departs 
and a new one arrives every .s days. When demand by each 
SSB(N) is assumed to be uniform, stock levels aboard the tender 
are depleted xt every s days. Theoretically, for the case 
where there is no re-supply, the tender is completely depleted of 


Spare parts at the end of Y days. 


a | 


Tender Support Cycle 
Pig.) 
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i, Routine Re-supply 
Consider the first factor influencing re-supply mentioned 
previously. The average percentage of the optimum amount (z: ) 
which the tender maintains will be defined as the effective inventory 
level. This level is a function of demand, re-order policy, and the 
lead time associated with the transportation of the item. Consider 
the case where demand for all items is uniform over a ninety-day 
period. Fig (6) is a plot of effective inventory levels at the tender 
versus lead time required for delivery as a function of re-order policy. 
The graph can be utilized in several different ways. Ifa re- 
order policy is chosen the graph indicates the inventory level that 
can be maintained for various items as a function of lead time. If 
lead times are fixed and a certain inventory level is desired to be 
maintained then the re-order policy can be obtained from the graph. 
It should be noted that due to the capacity of the ship and the 
lead time required for delivery of an item the tender will never 
average an amount equal to the load list quantity for that item. 
Several examples will serve to illustrate the use of Fig (6). 
(a) Consider the case where an effective inventory level 
of 90 per cent is desired for an item. The graph 
indicates that this level can be maintained by re- 
ordering when the stock level drops to 90% provided 


the lead time is 10 days. Another possibility is to 
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re-order when stock levels drop to 95% if lead time 
is 16 days. The re-order quantity in each case is the 
demand to date plus expected demand over the lead 
time. 

(b) If the lead time for delivery of an item is fixed then 
the re-order policy determines the effectiveness levels 
maintained. With a lead time of 12 days a 60% policy 
results in an effective level of 76%, an 80% policy 
results ina level of 84.5%, a 90% policy in an 89% 
level and a 95% re-order policy in a level of 93%. 

If the lead time is 30 days, which is analogous to that of surface 
re-supply by AK prehen the maximum possible effectiveness level is 
83% for the example eons 

It is obvious that high inventory levels can be maintained 
effectively only for the case when the lead times are short. This 
illustrates the requirement for air transportation of spare parts as 
opposed to surface re-supply. 

In addition to the effect on effective inventory levels maintained, 
air re-supply has the advantage of reducing the confusion factor at 
the tender when supplies are received. Rather than being saturated 
with parts delivered by an AK all at one time the tender is able to 
maintain better control over spare parts storage. This in turn tends 


to reduce ''loss'' of items due to misplacement. 


ait 





ES aww: 
ame emai a 
i 

—_—_— —— 

—>  —<—— A -  —d 





_ rr 
ee ee ee ee ee 


=_— ©@ eee" 
——— = i SS es = 9 


VS 





Re-order Policies 


OWT, peo 


IS 8F Of L2 


Gb zp 6€ 


96 Ef 





ADIJOq aepz0-a¥ so uotjoun J 
ese 
OUT} PEST “SA ST2AZT ssouaAaraTTY 






%09 


%08 
%06 
%S6 


ve 


Ud 


—-— SO ns wn en ew ewe eo 


ST 


cl 


\%) 


—=— PSS SS — i = 


6 


O€ 


OF 


OS 


D9 


OL 


08 


OOT 


%) S[PAST SSOUsATIITIGT 


( 


Fig. 6 


28 





ay, Determination of resupply method for NIS/Nonstandard items, 
Consider a spare part that is required by an SSB(N) and that is 
not in stock (NIS) on the tender. If it is assumed that a spare part 
NIS on the tender is always available at the stock point then the 
mode of re-supply transportation (air or surface) can be determined 


by using a criterion developed as follows: 


Let, Dax = Julian date of departure of AK from stock point. 
“VaAK = Days AK spends in transit to tender. 
S} = Days remaining until )“* SSB(N) starts patrol, 
ee Julian day that an SSB(N) generates demand 
) J 8 

for an NIS spare part. 

er = Days required for loading, unloading and 
installation of spare part. 

LT = Days required by spare part to travel by air 
from stock point to tender (MATS, commercial, 
etc), 


If R, +S; 2 Daxw+lan+C  . re-supply the spare part by AK; 
ff Ryo+ Sj < Daw +A +C >» re-supply the spare part by air, 
subject to the constraint Ry =- D5 2 Aa 

If Ry 34 on , then it is impossible to re-supply the spare 


part prior to the SSB(N)'s scheduled departure. 
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oe Conclusions 

1. The number of stockouts will be reduced if stocks are 
maintained in accordance with the model. The model maximizes 
the average number of demands filled over some length of time. 
This is equivalent to minimizing the number of stockouts over the 
same length of time. 

2. <As the spare parts stowage capacity of the tender is in- 
creased, the total worth (military essentiality) of the inventory is 
increased, but at a decreasing rate. For example, if a small 
tender and a large tender each have their cubic capacity increased 
by the same constant amount, the small tender has a larger per- 
centage gain in the military essentiality of cargo than the large 
tender. 

3. Items with high essentiality to cube ratios are generally 
carried at high stock levels. For those items with nearly the same 
pattern of demand (=) , the model maintains high inventories of 
high essentiality , low cube items, and low inventories at low 
essentiality , high cube items. 

4. When items must be deleted from inventory because of 
space limitations, those with the lowest essentiality to cube ratio 
are the first to be eliminated. 

5. The spare parts stowage capacity of the tender is the 


limiting factor in determining the length of the period over which 
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demands filled should be maximized. The load list should be designed 
to maximize demands filled over a given length of time. 

6. Re-supply by air is required to maintain high inventory 
levels of essential spare parts aboard the tender. In order to main- 
tain essential spare parts inventories close to the optimum load list 
level at all times, itis waaear to have short re-supply lead times, 
or a large tender capacity. Since tender capacity is limited, high 
inventories can only be achieved by utilizing the short lead times 
associated with air re-supply. 

Ue Re-supply ships cannot meet demands for non-standard 
items. The long lead time associated with ship re-supply precludes 
the possibility of meeting demands for non-standard items as they 
occur. For example, it is a rare occurrence when a submarine can 
execute an order for a non-standard item and have it delivered by 
ship re-supply before going on patrol again. 

8. The iterative nature of the model and the large number 


of items in a real life inventory make the use of a computer mandatory. 
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APPENDIX I 


Derivation of Optimal Tender Load List Model 


(1) ex. i = . ~S Fa) Axe, 
(2) bo C2) = s. (x2 VOX ¢ : 


Ze 
(3) SiCZY) = S_ KL dAXi ~ Zi im SNAKE. 
(4) ECxcj=s Seleor be C2). 

Essentiality of a spare part is an important factor in deciding 
on a load list which will tutes the average number of demands 
filled by the tender. On the basis of past usage data the model 
should indicate the quantity of each item to be stocked such that the 
total essentiality function M is a maximum subject to the tender 
stowage limitations. 

The function ay Xi 52 ) will be introduced to represent a 
utility valve associated with the ability of the tender to meet demand 


for an item. “Wile usta 2: 


a (Xi, Zi)= \ ; ae when ALCRz:e 
L= \,3 oeera Nl 

It should be apparent that if demand is less than what is 
stocked the tender is able to provide 100% service. On the other 


hand if demand exceeds the stock level then the tender does as best 


as it can by providing an amount Z<. 
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In accordance with the measure of effectiveness mentioned 
previously the objective is to maximize the total utility value as 
defined above 

wi 
Max, M(,,2,...2y)= BLZ gee zal 
subject to the constraint 
amen Se | Soe 
Nae 

ELZ SOK we. VA= Ne ves ste Gye Oc core VRE e a Idee cg 
where {(x,,x,--...-%s) is the joint probability function of 
Bes) 2s > es 
By assumption (7) this joint probablity density function can be 
written as: 

COX Kaees- Xe) = LUXOR Fn) 
i) 
[By definition J Siixddx p =i J 
As a result Ele Se (x: ,2Z2) | can be expressed as follows: 
ed oo 
=e WOKE zy, \ _ z a Oi uaF, Xi \dx a. ees Se Gra Ky Za) Fu lKnlaxy ; 

Let M: sila) = EL (x20) 


iM Po 
eee YN sidai\= § lk ZVEON AX: 


Maximize M\ (258es.. Zn) = = Wie CZ 
Subject to: = CZ Ci 
Xu 
The problem presented is one of maximizing a function of wu 


variables where the variables are related by a single constraint 
equation. The problem could possibly be solved by eliminating some 


of the variables using the constraint equation. Eventually the problem 
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could probably be reduced to an ordinary maximization problem. 
This procedure is not always feasible and therefore the technique 
of Lagrangian multipliers will be utilized. 

To find the points of we -“" oan eee se ee 2 where the 
variables 2,,2,---.-2n are related by the equation Zc2i-C= © 


the function Lis defined as follows: 


N 
= > Mi sila,y—-AT 2-6. 20— ern 


LEA 


where A isa Lagrangian multiplier 


>= 2 = ye (Misia )= ATE zi Ca) 


ae 


= a Lhiscavil— Aacz 
day V=2\,2,....N 


As a result,the following set of(n + l)equations in(n + 1) 


variables 2,,2,--- 2m and A 1s obtained. 


See Poy Gy 2 = 
3 (az.-Cys © 
ou 
3 


& 





Gcz ee ey =O 
20s 


By rearranging the terms the above equations can be put in the 


form 


i sdziy\ = Act 






aa Qe 
(6) SiGe \ \ 
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— equation (3) 
(6') sl &x XK FLO Vax + Zr : x dx. | = A( = == 
Lebnitz's rule is utilized to evaluate the left hand side of this 


equation since the units of integration depend on the parameter Z;: . 


Leibnitz's tr) 


“Rule ° Se) eC ecn = WU be) oa lb (4) == hla, tal) 


Yr) 
dhilx, 4) a 
7 cGy ot e 


eS KE OC NAX: +2 ae wax. | 
ae £C2)- — Recon Bo= Zo Qa f ak Ow VAXKL 
= t SLO) dx 
ae Ze ox 
Since + SG Oax = Ao = 5 Fay axe oe LOKI AX 
CS Ta 
Then J,, BOOdxL = AH J SK Ndxe 
By substituting and rearranging terms equation (6') reduces 
to the following - 
A— J Sx NOx. = AG) 
Ze 
es — ae Ge 
on . Si Gaadot = = \ Goad (5) 


Equation (5), subject to the constraint equation, provides the 
method where-by the set Z can be determined. The two equations 
determine a unique set of wales provided that for all = t (xK\>O 
for all values of Ac 

The distribution function , me) , of a numerical valued 


random phenomenon is’ defined as having as its value, at any real 
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number ZW , the probability that an observed value of the random 
phenomenon will be less than or equalto the number 2 . 

) / 

| (F(z) alee = 27) } -Since the left side of (5) is a distribution function 


the following condition applies: 


Off FuDax S41 , 





In order to solve for Zi in equation (5) we must assume a 
value for A .It would materially decrease the iterative require- 
ments if the possible values of could be confined to a portion 
of the real line. 

Zu 
Since O¢< i $. aj axn <4 
Th a le = 
en Ofi-dACS) 
—1<ArA(& M. 
Lert) => Me 2dr 
= Cis 
A~—rAC&) SA 
—A(&)40 = A20 
Mc 
: mM. 

These limits apply for all \ and therefore X must fall in 

the range. wor(Ss j >A2O in order for all items to be 
Cr 


included in the load list. 
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APPENDIX II 


Derivation of S:(2:) for the case where demand follows the normal 





distribution: 
ie oo 
Sez oe " KODAK + Zi . FX ae 
where ECx2) = tee. e7 -~+ Sys 


f, 0 Sx =L=S™ Fudd: +S. Fx) dx: 
ae 4. Lx JAX = 1 ae ree 
but Ne Siu ides = LA <\ 


- FO Nax. =f1- 11 -A( EY 
“~ Si C22) = f- Ke Foleydxt ~ Ze Ci C- ASS YY] 


In order to solve for S;¢z.) an explicit expression must be 


Ze 
found for a KL FiCK ANAK | The procedure is as follows: 








" + iy oe 
a1 § xX fo Vax: =f. oe Ax: 
+ (S34) > 
By taking the derivative of es c = =) the 
following expression is obtained: (haem) 
- > (%4-n2\* - (> 
ad LL 2. -) ee. se Xo ee 
Sime Var ee mee 
— X- Me) > 
To 
Vv Xr Me L. (X,- KM, 
ae Oe : Lee ey" j= aS ri) 
— ee ea 
AXIT G = 4 


Integrating both sides.of the equation between the limits of © 


e 





and “=, 
. A (Xo-™ -\ + an XN eS Ze “wc ae 
i aX; Var ~~ Waar SD 
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Rearranging terms: 
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Note that the expression on the left is the same as equation 
II-lexcept for a factor Sj . In order to obtain the same ex- 


pression multiply both sides by 9% which results in 
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Combining the above derived relationships an explicit equation for 


S72.) results: 
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APPENDIX III 

A computer program based on the flow diagram illustrated in 
Figure 7 was written in the Neliac language and compiled on the 
United States Naval Postgraduate School CDC 1604 digital computer. 
The Neliac statements are found in this appendix following Figure 7. ; 
Outputs for several choices of parameter values confirmed the re- 
sults shown inS ection 5, 

As indicated by the computer flow diagram in Figure 7 the 
first step in the operation of the program is to carry out the routine 
input procedures. 

The next step is the computation of the value of lambda equal 
to the smallest fraction, cs . Denote this minimum ratio as 


oF 


This value of lambda is used to compute all values 





of AREA= (1-H ) fori=1,...N. Each 2i is computed 


from: WG — eer F (EE)- AREA 
UrsA = DE eS 
\ awe 


Te exe (3S 


which provides an approximation of the normal distribution function. 


In the expression above: ef (=) 7 1 
where =\+0O. VAWAS2\ (= ez) + O.O88G4HO2'| (BY 
«06 0, 


+ 0.027 43349 a) — 0. COO39446 (2) 
+ 0. 003828915" (my 


IDetails of the Neliac language may be obtained from Prof. 
R.M. Thatcher, Operations Research Department, U.S. Naval 
Postgraduate School, Monterey, California. 
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The product C.c@i is computed, stored and summed. The result, 
a ZC is stored in the cell named SUMCZ. 

This sum, 22. » 18 then compared with the contents of 
Big C, the specified total cube available for the storage of spare 
parts. If SUMCZ > Big C, then the cube of the load list, for the 
given lambda, is greater than the capacity allowed for storage. When 
this occurs the item with the lowest essentiality to cube ratio is dis- 
carded, anew A = mints ) is computed from the remaining items 
and is used to determine a new SUMCZ. This new value of SUMCZ 
is compared with Big C. If SUMCZ remains larger than Big C the 
item with the lowest essentiality to cube ratio is discarded. This 
procedure is repeated until SUMCZ becomes equal to or less than 
the cube available, Big C. Whena new load list is found that meets 
the cube constraint then the essentiality of the load list is computed 
and a print-out is obtained for Ziy s@i5,A, Zz, and = M.S). 
A table of the items discarded is also provided. 

If the cube of the load list computed for the initial value of 
lambda ( A= wore me ) is less than the capacity allowed for 
storage, then a print-out is obtained for Zea A, SOR and &Mcs(z:). 
In order to find the load list that most nearly meets the cube constraint 


it is necessary to evaluate 2c.zZi for several assumed values 


of r 
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The number of values of A to be assumed is specified in the 
program by the parameter K. If, for example, K=20 then the 
variables <v, S(2.) 22 and &Mis(z.) will be evaluated for 

| _ My 
twenty equally spaced values of X between min — and zero. 
In this case twenty print-outs will be provided. The optimum load list 


is taken as the load list that most nearly satisfies the cube constraint. 
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PARAMETERS 





Min = —> MIN 
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“DAS TRABATTIION 


COMPUTES 





TSB. SOMCH] 
—e SUMCZH 


sar A 











a 

For L=NG to N-1 Compore ™ SG; 
ITEM 14 et + SomMmsz 
(inciwoEs Mi, —F SoOMmMS Ss 

Se ymi, TE ) 

N-4-—N 
DISCARD TABLE 
Fig a Computer Flow Diagram 


42 





“PRINT (TEM 


DISCARD 
TAS\& 





FANIS® 


— MIN 
LAMBDA a 





FINISH 


‘Fig 7. . Computer Flow Diagram 
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Neliac Statements for Computer Program 


Big C= 1084 1*¥O 


NWN = \6, 
MX Gere 100.0 *¥ OD, 190.0 * 0, 
eerie SO.O * O, 
\OO-O ¥O, 100.0 * O, 
FO.0 *¥O, SO-O# DO) 
\OO.O ®O, 190-O%O," 
SO-O *O), SOOXO, 
\OO8-O KO, IOO.0 KO, 
50.0 *O , 0050.0 *O, 
A 
SIGMACG) = \2.0KO 5 OO KOS 
Oe, 2G rn, 
\A.O *O y BO *O, 
R2OKD, 3%O *O, 
LO KO, BOK OY 
Weoe~ © BORO, 
20 KO, 30 ®O 
12.0 ¥O , CO3.0 #0, 
WG) = AOKO, AO KO, 
A.O KO} AO KO, 
AQ KO ,y AO XO, 
AO KO yg AO *® Oy 
ok >). +O*X OY 
A2AO*OD 2ZAO KO 
2zO KO, BOKOY 
ROXG, COOZ2.0O *XO*= 
CUS) = WOKO , WOKO, 
WO FOL WOK ©,*7 
2O KO 5 3.0 * OY 
AO KO , WOKS, 
(ae 5 WOO" 
oO KO , 3.0 *O, 
BO RO 5 OCOR.OXKO, 


ITEM NRCG) = 4,2, 3,4,5,6,17, 8,9 


> 
WO, 11,125 15, VA, Sey 4S 
MIN = © x*O 


KK = 20.0 *O 


ae O00 , 
LAMBDA = COOO.0Cccd *O : 
SOMCe = OOOO. COSO * CO, 


SOMNISZ. = QOOCCO.50O05 * or, 
AREA Ube) = Ox © 


Ze) = Qcocce. Sac xO, 
S26) = OGO000 . Coo * O, 


NN. Indicates a blank card 


9 Indicates a card of all nines 
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Reco. = Gerke... 


RSPAS = O50. OG ta Fe, 
Nears = — = *¥ O 
Pee DS Oe. N = OSoreoce -& © 
NGb2exr = O #* © 
GxPi = ©x CQ, 
2xXr 2 = ero. 
allio = a= 
TI NIipeaiie= alos ZassSsonz x Oj eee, 
NN 
MELAG , me , Mie mee, 
> 


: SELAG SVM, DISCARM C\L)\ = O00 5 
ea Mu en auc) 
Sok = 


3 
= PVT AT COONTININ OSD 2 5 
A 
3 SY, Ca ae 


; i Pe PROGRAM | TAMICD \ HARNEY 974, , i 


PRIRYT VEE Rte 
AN 
Om NG > FLAG — TTLAG —e MFLAG —S N, 
SQRTE CVYRbPT ; SQATZPL) , 
i 
FIND Mins “WoC : 
Og of (ose. 
REPEAT: 
TF Tem NR Tma = NOU. M+ Lem, MeN: 
; ANUS tA. REPERT, ° 35 
Wier / ¢ Om = LAmena 
éN 


OC —— ae MIN FLAG, 5 
FER K = MEV) Ny -— A | 
{EF Wem ne OKT] BOUL: DISCARD Wb. 4 
WORT / CORY —e Wa Ae < CAMBDA : 
MID C — LAMSEDA >. 


Ke Min FLAG 5; 
ASCARD Wid i 
A, 


Tr MELAG = Mi LAMBDA —— Mrnl, 1 — MEAG 
m, 
COMPUTE. AREA + 


On > SemMce , 
FOR, W= QCA) SH- 1 


{ XF Trem wea Cid = SUL 


Bohrer 

) 

V.O— LAMEDA & (CUT) [ws ) > reea CE), 
Been CL \=S 


2 
~~» 


On 3-87 — > TT COn™. ;: 
TF aAnen LEA 0.5 1: 0.5 —Arca EX 


4 

—> ARTAMS 

ARZAMS &— 0.000024 + OT; CH, 3 ? 
NORMEX CAREAMS 3¢ a M9 

necmt PEARS OPS RENE 

R NORMEX CAREAMS 2-7 13 


CibnyT 
SD 4 Aer Ogee 
» i << ocx ist = Ix <4 LAMBOA \ = \ 
> LAMBIVA < WL) AacaA | = 


\ > AREER LE) <. 


45 





y 


2 7 2 


UW AREAMS | 
oa Oke VOT l= \ oy > 


Mx CLA ~ SIGMA CT) KT — =2CTI ; 
eae K CUS \ + SOMCe —— SUMCeH , 
BMvr és s 5 


» L PRINT << Sumc#® l=) > suomcz >3, 
TF JFLAG =O . TEST, CPMePoTE 8z. 


VEST. SOmMcH# < BAGC :. CHOmeclesz. 5 


4 —» FLAG, 

NG~L —e NG 

MIN FLAG —e Lo = ITEM DiscanD Enid , 
N+. —N, 

NOU —~ TEM NR EWU), 

FINnDo MIN VIDC. 


CHMPOTE SE * <2? 
, Cec bun & UG’ste >>, 5, 
Or SOMWSS , 


ON: FER T= OOUNN—L LIF TEM NRE) = NOLL: SKIP; 
SIGMATETA /SQRTZP©R — Sieh2zer, 


MxXTTA / SIGMA CI] —e KKK —™ XK, 
NEGETS * X—> Ky, Expr (XX, EXP), 


(ZCTIA— MXCTIA/S\GMATI] — VX YN, 
NEGrr S «VY — VY, Sxer (VX 3; =Xpz2), 


MXCUEN & ARCA CE] + ZCE] * G.0 —Aren GI) 
+ SAGH2PL *® EXPLA — ANGOzST * EXP2—K—S= ,y 


WEE) *« SH ~*~ SOMWSSZ FP SOME, 
ce A —— TL 
[ aw << (NNN “eG Me see > 3, SKIES, 


DASCARD TABLE - 

3 L Sena ay 9-9 3 

5 Y << rrem \ discaRD | TABLE. 77 959 $ 
5 EK WIT OM SN 

OY CG UNAM UY SIGMA 


FOR T= @ CA NG — AL 
{ rem Msenayws CFl— JS, 
TEM Dscaew TKI + 4 TEM Dscaawy Cr) 
{ PawT < \ TEM Discard CLA Ww CJD \\V 
— Casa Mx CI] WW HGMACTI > $, oe 
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wor. 
a. 











2S Uew~ [oie = . 
, street ee 1 © 
623. eae ™ 
- » Jae, = ow, ae eee, 
OS eee, te .) 
(aa o — Pane ee “sa 4 @ 
a ae J* 9/0 oo pa a hw 
(fe gee 6 “eet «- ee & fe . 
eee — “he! @ « - -e 
a a - 
i. =? ( Ss « 6 eat —— | es 
iia: 
'oe ii “ke 
o —————— 1, [= 
\— Lao = 
= -& ; . 
,™ 
J “ais ios * <= 
:¢ SS ——eeeee alll 
ii _—_ * 1 => =—__ ee 
a - a"; _ - 
‘ee? = | 





A 
3 


Op? 


A 


FIniseH <% 


=) 
JJ = CO 


> PRT PARAM : 


L 
, Lo MET YVUWUCUM WISGMA >, FY 
a >. 


aa <j — WWD WC CTT WW MX CTT HGMATTI 
> % : 
3 3 = 


5 

A\=O x*O, B\VFTOKO 
D \ SuTrarK, = O*% © 
N \ 
= 


OM = O*XGD, VIEW Am = ORO, 
LRoTZ = VMHANRZVESS XO 

RA\GPT2APE = ©. 39594228 * Oy 

CI RE * 


Ngdemex (AREA. TT) +2 | 
TF AER EO 1. OFT, EXETRX . 
Ir OF AREA WH AREA y o.S > SRAM? . 5 
Ve) = 
Clameucre » V/S8T2 —e A, ERAGAF (A, ERE), 
QS * ERE — ARCA — PF NOM, 
AN *khmmA.L, Om-A HA : 


REGTZAYPIO x A ——F DENDM , 
NOM /DENADAM —F>_DELTATK , 
TT DECTATK VG 
= aS). So te anne oe 


WELTATR, fT B&B 
O-—- ©. OocoSs <~H< O.00005 4 


SXETROK. CPMPTE | 
Glarge i. , 146 Tireeac_lwvuave | 


OE | AwEaA \ WW | 
mfemex | Fosenda 27 5, 


=e \Xrmx:$.. 


A 
2 


) 
ERRORE ON. ERE. P.) + 


L 6660 &. COZZEAMIS KF WG. COS SAWN) 
% WN 0. 02-74 B3B4Q, ) HW 1 O.08 564027 ) 
XK W+ANIQUABAL) *¥ Wr s.o ——P, 
Px P—~P, Ww x*« DP —~ Ps 
i oe 
4.0 —1.07 p —~ ERF,4.. 
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pow 


nn 


zee i, 
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